library(dplyr)
library(tidycensus)
library(ggplot2)
library(mapdeck)
library(shiny)
library(shinythemes)
library(shinyjs)
library(ggplot2)
library(maps)
library(plotly)
library(DT)
library(dplyr)
library(tigris)
library(tidyverse)
library(tidycensus)
library(readxl)
library(collapsibleTree)
library(shinycssloaders)
library(leaflet)
library(leaflet.extras)
library(rvest)
library(sf)
library(shinydashboard)
library(shinydashboardPlus)
library(tidygeocoder)
library(janitor)
library(dplyr)
library(ggplot2)
library(magrittr)
options(tigris_use_cache = TRUE)
library(ggwordcloud)
require(rgdal)
ind <- st_read(dsn = "/Users/nandinidas/Desktop/2022-DSPG-LivDiv-/data/GADM", "gadm36_IND_3", stringsAsFactors = TRUE)
## Reading layer `gadm36_IND_3' from data source 
##   `/Users/nandinidas/Desktop/2022-DSPG-LivDiv-/data/GADM' using driver `ESRI Shapefile'
## Simple feature collection with 2340 features and 16 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 68.18625 ymin: 6.754256 xmax: 97.41516 ymax: 35.50133
## Geodetic CRS:  WGS 84
sundarban <- subset(ind, NAME_2 %in% c('North 24 Parganas','South 24 Parganas'))
d.sundarban<-st_union(sundarban)
#summary(shape@data)

#shape$INTPTLAT <- as.numeric(shape$INTPTLAT)
#shape$INTPTLON <- as.numeric(shape$INTPTLON)

#par(mar=c(0,0,0,0))
#plot(d.sundarban, col="darkgreen", bg="white", lwd=0.25, border=0 )

#shp_df <- broom::tidy(shapemap, region = "")
#lapply(shp_df, class)
village_all <- st_read(dsn = "/Users/nandinidas/Downloads/fgd-village-selection-main/shapefiles", "Village, GP coordinates", stringsAsFactors = TRUE)
## Reading layer `Village, GP coordinates' from data source 
##   `/Users/nandinidas/Downloads/fgd-village-selection-main/shapefiles' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 1054 features and 11 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: 88.04613 ymin: 21.56433 xmax: 89.03389 ymax: 22.86735
## Geodetic CRS:  WGS 84
village <- subset(village_all, Village.Na %in% c("Amrabati","Beguakhali","Bijoynagar","Birajnagar","Haridaskati Samsernagar","Lakshmi Janardanpur","Pargumti","Purba Dwarokapur","Sagar","Shibpur"))

icons <- awesomeIcons(
  icon = 'ios-close',
  iconColor = 'black',
  library = 'ion',
  markerColor = "lightred"
)

leaflet(data = d.sundarban) %>%
  addTiles() %>%
  addPolygons(
    fillColor = "green",
    stroke=TRUE,
    weight = 1,
    smoothFactor = 0.2,
    opacity = 1.0,
    fillOpacity = 0.7,
    highlightOptions = highlightOptions(color = "white",
                                        weight = 2,
                                        bringToFront = TRUE)) %>%
   addAwesomeMarkers(~lon, ~lat, label = ~as.character(Village.Na), icon=icons, data=village)

#Birajnagar is right behind Bijoynagar

Pragumti is not there — I found Praghumti

Sagar is not there too — Gangasagar and Sagar Madhabpur